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Abstract. We discuss the properties of devices of small Josephson junctions in the light of the phase 
fluctuation theory and the energy band structure, which arises from the derealization of the phase variable. 
The theory is applied in the realization of a mesoscopic amplifier, the Bloch oscillating transistor. The device 
characteristics and comparison with theory and simulations are discussed. The current gain of the device 
in a stable operating mode has been measured to be as high as 30. Measurements on input impedance and 
the power gain show that the BOT is an amplifier designed for middle-range impedances, ranging from 
100 kfi - 10 MO. 

PACS. 73.23.Hk 74.50.+r 73.23.-b 



, ^-v 1 Introduction 



The conjugate nature of charge and phase gives rise to 
a rich array of physical phenomena, which have been in- 
tensively studied both theoretically and experimentally. 
The quantum nature of the phase variable was shown in 
macroscopic tunnelling experiments |1] and its conjugate 
relationship to the charge has been shown in many conse- 
quent studies [SJ. 

One of the consequences of the charge-phase conju- 
gate relationship is the Coulomb blockade of Cooper pairs 
which arises in very small Josephson junctions |3J. The 
system is described in phase or charge space, depending 
on which variable is a good quantum number and, con- 
sequently, determines the eigenstates of the quantum sys- 
tem. Associated with the two variables are two competing 
energy scales, the charging energy Eq = e 2 /2C and the 
Josephson coupling energy Ej. The Hamiltonian for the 
small Josephson junction has a periodic potential, hence, 
Bloch states and a corresponding energy band structure 
can be derived — analogously with the conduction elec- 
trons in solid state physics |4l5j . 

Incoherent tunnelling, the interaction of tunnelling elec- 
trons or Cooper pairs with the electromagnetic environ- 
ment, is a strong effect observed in small tunnel junctions 
both in the normal and superconducting states |6l7j . Small 
tunnel junctions can be used as sensitive detectors of envi- 
ronmental impedances and noise sources but, this sensitiv- 
ity also means that attention to control of the electromag- 
netic environment is essential when designing mesoscopic 
devices. 



The Bloch Oscillating transistor is a mesoscopic de- 
vice which is based on the dynamics of the Bloch bands 
in a voltage biased Josephson junction (JJ) in a resistive 
environment. The main operating principle was demon- 
strated in references [81911011 T] . The conclusions that can 
be drawn from these articles is that the BOT shows con- 
siderable current and power gain (for both cases ~ 30). 
Additionally, the noise temperature at the optimum oper- 
ation point was observed to be as low as 0.4 K, although, 
the theoretical and simulated noise temperature has been 
shown to be 0.1 K or even below. The BOT has an input 
impedance, which can fairly easily be tailored by fabri- 
cation to be in the range 100 kSl — 10 Mf2. Thus, the 
BOT is a device suitable for use at mid-range impedances 
and it complements the group of basic mesoscopic devices: 
the single-electron transistor (SET) - a device for high- 
impedance applications and the SQUID, a low impedance 
device. 

At this point, we would like to emphasis that one 
should not confuse the BOT device with the Bloch transis- 
tor [12], which is a similar but, fundamentally, a different 
kind of J J device where the current of a single Cooper pair 
transistor (SCPT) is modulated by a gate voltage. 

This paper will briefly review the physics that gives rise 
to the voltage-current relation of the Josephson junction in 
light of both the Bloch band model and incoherent Cooper 
pair tunnelling. After the theory, a few observations on 
the "Bloch nose" - a consequence of Bloch Oscillations 
- will be presented. The last part of the paper will take 
on a more device oriented approach with the discussion 
of the experiments and the question how well the findings 



2 R. Lindell, L. Korhonen, A. Puska and P. Hakonen: Incoherent Cooper pair tunnelling and energy band dynamics 



can be reproduced by theoretical simulations on the Bloch 
Oscillating Transistor. 



2 Theory 

2.1 Theory of Bloch states 

In mesoscopic systems the phase and charge comprise a 
conjugate pair of quantum variables, analogously with the 
more familiar space and momentum coordinates. The vari- 
ables thus satisfy the commutation relation 

[4>,Q]=2ie. (1) 

The total classical Lagrangian for a single, isolated Joseph- 
son junction can be written as 

£= < ^+E J coscj), (2) 

which consists of the charging energy due to the capaci- 
tance C and the Josephson coupling energy as the poten- 
tial. From the commutation relation ([1]) we can immedi- 
ately write the quantum mechanical Hamiltonian [5] as 

H = -fi C g^,M 2 - Ej cos 4> - —I<t> + H env + H mt , (3) 

where also the interaction between the driving current and 
the phase variable oc — lip, the junction environment en- 
ergy H env and the coupling with the environment Hint 
have been included. Depending on the theoretical frame- 
work used, some, or none of the last three terms are in- 
cluded. 

To set the initial stage, we consider only the first two 
terms in the Hamiltonian, in the case when Eq 3> Ej and 
the charge is a good quantum number, thus leading to 
Coulomb blockade of Cooper pairs and a complete der- 
ealization of the phase. Equation ^ then takes the form 
of the familiar Mathicu equation with the well-known so- 
lutions of the form &%((p) = e ivq ^ 2e u n ((p), where u n (ip) 
is a 27r-periodic function and the wave functions are in- 
dexed according to band number n and quasicharge q |13j . 
More precisely, the Mathieu equation has also 4ir periodic 
solutions, which correspond to the the case of quasipar- 
ticle tunnelling [14j . The resulting energy band structure 
is illustrated in Fig. O Verification of the existence of the 
energy bands has been carried out by different methods 
|15 | 16 | 17| . 

For the opposite limit, Ej ^ Ec, we should index the 
system eigenstates according to the phase variable and 
the charge would be completely delocalized. This situa- 
tion corresponds to the classical superconducting state of 
the Josephson junction, which can also be described in 
the "tilted washboard" picture when including the third 
term, cx Icf> in the Hamiltonian ([3]) . This model can be eas- 
ily extended to the well known resistively and capacitively 
shunted junction (RCSJ) model, where one includes a dis- 
sipative resistor in the circuit. But, we will not discuss or 



apply the RCSJ model in this paper as it is mostly used 
in the case of large Josephson junctions in the context of 
macroscopic quantum tunnelling. 

Next, we consider the Josephson junction in the band 
picture, when the phase is delocalized and the state of the 
system can be described by its quasicharge q. For a more 
realistic situation, we also need to consider the junction's 
coupling to its environment. The detailed and involved 
calculations for the cases with ohmic and quasiparticle 
dissipation can be found in Refs. [5114) . The physics is 
rich in detail and includes phenomena such as the zero- 
dimensional quantum phase transition when the system 
goes from a superconducting R < Rq to an insulating 
state R> R Q [T5]. 

The instantaneous voltage over the junction is given 

by 

v^-EM. ( 4 ) 

Aq 

We consider the system to be either current or voltage 
biased, meaning that the environmental resistor is situ- 
ated in parallel or in series with the junction (see Fig. [I]). 
With a steady current / flowing through the junction the 
quasicharge evolves according to 




Fig. 1. Schematic view of the series or voltage bias (left) and 
the parallel or current bias (right) configurations as defined in 
the text. 

In the voltage biased case, the voltage over the total 
system consists of the voltage drop over the series resistor 
and the junction voltage 

V x = IR + V. (6) 

In the current biased case, where the resistor is parallel 
with the junction, we have the relation for the total cur- 
rent 

I x = V/R + /. (7) 

Hence, one can switch between these two cases if one de- 
fines I x = V x /R 

In practice, we fabricate the resistor in series with the 
junction, and thus consider next the former case of voltage 
bias. If the resistor R is large (compared to the CB resis- 
tance of the junction), we can still think of the junction 
itself to be current biased. Hence, if the driving current 
is low enough, Aq/At <C e8Ei/H, where 8E\ is the gap 
between the first and second band, the quasicharge q is 
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increased adiabatically and the system stays in the first 
band. We are then in the regime of Bloch oscillations; the 
voltage over the junctions oscillates and Cooper pairs are 
tunnelling at the borders of the Brillouin zone, or as for the 
definitions here (see Fig. [J), at q = ±e. Consequently, the 
current through the junction is coherent and the voltage 
and charge over the junction oscillate with the frequency 



Sb = I/2e. 



(8) 



The theoretical I x — V characteristics for an external 
current bias thus first shows an increase of the junction 
voltage with increasing current but at the onset of Bloch 
oscillations, the voltage decreases. If the current I is not 
adiabatically small, we can have Zener tunnelling between 
adjacent energy bands. The tunnelling is vertical, i.e the 
quasicharge does not change. The probability of Zener 
tunnelling between bands n — 1 and n when Ec 3> Ej 
is given by 



P n.n-1 = CX P - o 



7T 5El 



nE c hi 



= exp -- ^ , (9) 



where SE n — E n — E n —\ and Iz is the Zener break down 
current [TO] , 

The downward transitions can take place through sev- 
eral processes. The cases we need to consider are transi- 
tions due to quasiparticle tunnelling and due to charge 
fluctuations caused by the environment. The quasiparti- 
cle transitions in the JJ couple q states that differ by one 
e. In the BOT, the transitions are primarily driven by 
a current bias through a second junction, explicitly de- 
signed for returning the system to the lowest band. The 
external environment gives rise to e.g. current fluctuations 
that couple linearly to the phase variable. These can cause 
both upwards and downwards transitions. The strength of 
the fluctuations is given by the size of the impedance: the 
larger the impedance the smaller are the current fluctu- 
ations and the transitions rates. As we will see later on, 
the successful operation of the BOT requires one to con- 
trol both the upwards and downwards transition rates, 
or rather their relative strength. Later, when discussing 
a very simple BOT model, we will make use of the Zener 
transition rates and transitions due to charge fluctuations, 
both derived in Ref. [Ml. 



2.2 Incoherent tunneling and phase fluctuation theory 

It is well known that the electromagnetic environment of 
tunneling junctions affects the tunneling process by allow- 
ing exchange of energy between the two systems formed by 
the tunnel junction and the external circuit |6|20|21|2"2] . 
The influence of the external circuit can be taken into 
account, for example, within the so called P(_B)-theory, 
which is a perturbation theory assuming weak tunneling. 
The external environment is thought of as consisting of 
lumped element electric components which give rise to 
a classical impedance Z(ui). The impedance is quantized 
using the Caldeira-Leggett model, which consists of an in- 
finite number of LC-oscillators. Hence, in this limit, one 



can model dissipative quantum mechanics. The energy fed 
into the bath never returns to the tunnel junction system. 

A perturbative treatment of the Josephson coupling 
term gives rise to a simple looking result for incoherent 
Cooper pair tunneling |2QI22j . The forward tunneling rate 
is directly proportional to the probability of energy ex- 
change with the external environment: 



r(V) = —E}P(2eV), 



(10) 



and the backward tunneling rate is r (V) = T(—V), thus 
leading to the total current 



I(V) = 2e (r(V) - r(V)) 

= (P(2eV) - P(-2eV)) . 

The function P(E) can be written as 
1 



(11) 



P(E) = 



2irh 



dt exp 



J(t) 



-Et 



(12) 



which is a Fourier transform of the exponential of the 
phase-phase correlation function 



j(t) = <[*>(*) - v(Q)] v(o)) 



(13) 



The phase-phase correlation function is determined by the 
fluctuations due to the environment and it can be related 
to the environmental impedance with the fluctuation-dissipation 
theorem. The result is that J(t) can be found by 



J(t) 



dLuReZ t (uj) 
u R Q 



{coth(/3fiw/2)[cos(o;i) - 1] 
— i sin(ajt)}, 



where Rq — h/4e 2 , and 



)C+1/R 



(14) 



(15) 



is the impedance seen by the junction. 

We also need to consider quasiparticle tunneling, es- 
pecially in the BOT circuit, where we also have an NIS 
junction. We will not, however, take into account Cooper 
pair transfer by Andreev reflection [23] as it is not essen- 
tial in the experiment. For quasiparticles the perturbative 
tunneling rate assumes the form 



r„ 



i 



e 2 R T 



dEdE' 



,Afi(g) N 2 (E' + eV) 
JVi(O) JV 2 (0) 



(16) 



xf(E)[l-f(E' + eV)]P(E-E'). 



Here, the density of states on the two sides of the junc- 
tion depend on whether the lead, i = 1 or 2, is normal or 
superconducting. In the superconducting region, we have 
Ni(E)/Ni(0) = \E\/VE 2 - A 2 , for \E\ > A and zero oth- 
erwise, and in the normal region Ni(E)/Ni(0) = 1. 
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2.3 BOT conceptual view 

The circuit schematics and a conceptual view of the BOT 
operation cycle are shown in Fig. [2J The basic circuit ele- 
ments are the Josephson junction, or SQUID, at the emit- 
ter, with a total normal state tunnel resistance of Rjj, 
the single tunnel junction at the base with the normal 
state resistance Rn, and the collector resistance Rc- The 
BOT base is usually current biased with a large resistance 
i?s at room temperature, but the large stray capacitance 
Cb results in an effective voltage bias. A requirement for 
the successful BOT operation is that the charging energy 
Eq — e 2 /(2C) is of the same order of magnitude as the 
Josephson energy Ej. For the theory presented in this 
paper to be valid, the condition Ej/Ec <C 1 has to be 
satisfied. 




Base 







sm< 1 







q/e 



Fig. 2. Circuit layout of the BOT and conceptual view of its 
operating cycle. The island is marked by O- The two band 
approximation of the BOT dynamics: Bloch oscillations (red 
arrows), Zener tunneling and Coulomb blockade (black arrows) 
and relaxation due to quasiparticle tunneling (black dashed 
arrow) . 



The theory describing the workings of the BOT is out- 
lined in articles [1018111] . The work by Hassel and Seppa, 
|llj . contains both an analytical approximation of the 
physics and the basic principles for a simulation which 
are also used to some extent in this paper to test the 
agreement between the experiment and simulation. 

The basic physical principle of the BOT relies on the 
existence of Bloch bands in the J J, which is embedded in 
a resistive environment Rc ^> Rq = h/Ae 2 . The BOT 
is voltage biased to a point on the lowest band where 
Bloch oscillations can start. The supercurrent thus flows 
due to coherent Bloch oscillations (see Fig. [2j in the low- 
est band until the flow is stopped by Zener tunneling into 
the second band. The JJ is Coulomb blockaded until it 
relaxes down to the first band, either intrinsically due 
to charge fluctuations caused by the environmental re- 
sistance or with a controlling quasiparticle current. The 
intrinsic relaxation is detrimental for BOT operation and 
thus the fluctuations should be kept low by requiring that 
Rc ^> Rq = h/Ae 2 . In practice, we need Rc > IOORq 
to be close to the ideal operation as outlined here. The 
control current is injected through the second junction, 
which in our case is a normal-insulator-superconductor 
(NIS) tunnel junction. The amplification mechanism can, 
in its simplest form, be said to arise from the number 



of Bloch oscillations triggered by one quasiparticle. The 
Coulomb blockade of Cooper pairs (CBCP) is thus a ne- 
cessity for the BOT to function. 

Our simulation of the BOT is based on the method 
given in Ref. llj. In the simulation it is assumed that 
the current flowing in the different parts of the BOT: the 
JJ, NIS junction and resistor can be treated separately 
and the dynamics of the island charge (Qj = ViCjj, 
where Cjj is the capacitance of the Josephson junction 
and Vi is the island voltage indicated in Fig. [2J is simu- 
lated and averaged over a large number of steps (typically 
10 - 100 million). The tunneling rates or, rather, the tun- 
neling probabilities in the JJ and NIS are calculated us- 
ing the P(i?)-theory [22] described earlier. Furthermore, 
we assume the same resistive environment for both junc- 
tions. The capacitance in Eq.[l5jis then C — Cjj + Cnis- 
The tunneling at each point in time is then determined 
by comparing a random number and the tunneling prob- 
ability of the different junctions, which depends on the 
voltage across them. Furthermore, the simulation actu- 
ally assumes that Ej / Ec <C 1 and that the energy bands 
can be approximated by parabolas: E(q) = q 2 /2C. In this 
case, the quasicharge is equal to the real island charge 
Qi, which is, according to the model in Fig. [5] also the 
charge over the Josephson junction. The equation for the 
simulated island charge is then given by 



dQi _ V c ~Vi 
dt R c 

I dt 



V dt /Q/V< 



V dt 



QPj 



(17) 



CP 



Here the collector voltage Vc takes the role of V x in the 
theory of Bloch oscillations considered earlier and Vi is 
the island voltage. The island charge is thus modified by 
four terms: the constant relaxation current through the 
collector resistor, the quasiparticle tunneling through the 
base junction (QPnis) and JJ (QPjj), and the Cooper 
pair tunneling in the J J. The theory thus excludes any 
quantum mechanical interactions between the Josephson 
and NIS junctions. This simplification has shown to be 
quite useful in determining the main operation princi- 
ples of the device. Although, one can argue that this sim- 
ple treatment of the island voltage as a time dependent 
variable while the tunneling rates are unaffected by is- 
land dynamics is not the correct way. As pointed out in 
Ref. [23] , the average fluctuations of the phase of the island 
will govern the tunneling in the Josephson junctions and, 
therefore, a more rigorous approach would account for the 
peculiar phase fluctuations from the tunneling quasipar- 
ticles. Nevertheless, time dependent P(£?)-theory should 
provide a reasonable starting point when the environmen- 
tal impedance is large [25] ■ However, how to include the 
Zener tunneling and the dynamics of the Bloch bands into 
the phase fluctuation model is still under investigation. 
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2.4 Alternative analytical theory 



where v — CVc/e, r = RcC and 



The physical principle of the BOT current gain can also 
be derived analytically from another viewpoint [9]. The 
average BOT emitter current can be thought of as the 
result of being in either of the two states: the Bloch os- 
cillation state with a time-averaged constant current and 
the blockaded state with zero current: 



(18) 



r Vc/Rc, n = i/r T 

1 o, n = i/(n+r B ). 

The amount of time the system spends in each state is 
given by the Zencr tunneling rate, If, the intrinsic relax- 
ation If , and the quasiparticle driven relaxation Tb ■ The 
base current, however, flows during the opposite times: 



o, r T = i/r T 
e r B , n = i/(r l + r B ) 



(19) 



From these equations we can simply derive the average 
emitter and base currents 



(Ik) = 



(Ib) 



Rc n + t l 
n 



(20) 



(21) 



From Eq. (f2Tj) we can solve for Tb and insert this into 
Eq. (|20p in order to find the emitter current 



(h 



v c r. 



V c 



1 



R c r T + r L eR c r T + r L 

We thus find the current gain 

d (I E ) V c 1 



(Ib) 



d (I B ) eR c I> + If 



(22) 



(23) 



The base current relaxes through the collector resistance 
and, therefore, the collector and emitter gains are related 

by 

0e = 0c - 1. (24) 

The current gains are defined with the minus sign for con- 
venience, because, both theoretically and empirically, with 
the sign convention used here the derivatives are negative. 

Next, we have to find the transitions rates If (Vp) and 
If (Vc) as a function of the collector voltage. These have 
been calculated by Zaikin and Golubev [35]. The Zener 
tunneling rate in a resistive environment, and with the 
assumption Ej <C Ec, is given by 



r T = -exp 



vz 



1 



{Sg 2 /e 2 ) 
(v-1) 2 



(25) 



and the down relaxation rate due to charge fluctuations is 
given by 



r, = 



: exp 



(«-l) 2 

2{5q 2 /e< 



-{v-lf 



(26) 



vz 



TT 2 R C 
8Rn 



Ej_ 

E c 



(27) 



The voltage vz is naturally related to the Zener break 
down current Iz = evz/{4=T). Equations (f2"5| and ([2Tj)) 
are given as function of the bias voltage v, although, the 
theory is strictly valid for the current bias case. One can 
switch from the current biased to voltage biased case by 
doing the transformation as discussed earlier by setting 
I = Vc/Rc- In our experiment, the relevant case is the 
voltage biased one. 

Closed forms for the charge fluctuations in a resistive 
environment can be found in the two limits of thermal and 
quantum fluctuations: 

r k B CT, k B T > a s E c 

^ 2 ln{uj c T), k B i C Q s -fc c , 

where a s = Rq/Rc- At the low temperature limit where 
quantum fluctuations dominate one needs to know the cut- 
off frequency u> c , which depends on the details of the ex- 
ternal impedance Z(u>). In practice, our experiments are 
done in a regime where ksT ~ a s Ec, where both effects 
are present. 

The two-state model presented here has also been used 
for the basic mechanism that generates the output noise 
of the BOT [11]. The noise of the device is found to be 
output governed and the input current noise is, therefore, 
inversely proportional to the current gain [27] . 



3 Fabrication and measurement 

The fabrication of the BOT is done using standard elec- 
tron beam lithography and 4-angle shadow evaporation. 
The order of the evaporation process was chromium (1), 
aluminum (2), oxidization (3), chromium (4), copper (5) 
and aluminium (6). Originally, the BOT device was thought 
to have an NIN junction as base junction, but the tech- 
nique of manufacturing both SIS and NIN junctions on the 
same sample has not yet been mastered, therefore, we used 
NIS junctions instead. In the first design, the NIS junction 
consisted of the aluminium-aluminumoxide-copper inter- 
face, but it turned out to be the most sensitive part of the 



BOT # 


Rn 


Rjj 


Rc 


Ej 


Ec 


Ec 


Ec 


1 

2 


90 

53 


24 
8.8 


188 
368 


0.28 
0.37 


0.93 
0.78 


0.3 
2.1 


0.03 
1.8 



Table 1. BOT parameters for two measured samples. Rn and 
Rjj are the normal state resistances of the NIS and JJ tunnel 
junctions, respectively. Resistances are given in units of kf2 
and energies in Kelvins. The tunability of the Ej/Ec ratio is 
indicated by the max and min values. 
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device and it often was already broken after bonding to 
sample holder. An additional, 7 nm thick chromium layer, 
seemed to protect the junction and improved the sample 
yield considerably. 

Instead of having a single JJ at the emitter we made 
two junctions, which thus form a SQUID loop. The SQUID 
configuration gives the possibility of tuning the effective 
Josephson coupling energy according to 



Ej 



Ej 



2EjEj cos(7T$/$ ), (28) 



where Ej and Ej are the Josephson energies for the two 
junctions, $ is the externally applied magnetic flux per- 
pendicular to the loop area and <&o = h/(2e) is the quan- 
tum of flux. The Ej tuning allows us to find the optimal 
operation point in terms of the current gain [8] as well 
as giving the possibility to compare the Ej dependence 
with the theories presented in this paper. In practice, the 
asymmetry of the SQUID junctions makes the tunability 
less than perfect as can be noted from Table 1. 

The BOT measurements were done on two different di- 
lution refrigerators: a plastic dilution refrigerator (PDR- 
50) from Nanoway and the other from Leiden Cryogenics 
(MNK- 126-500). Both had a similar base temperature of 
30 mK. The filtering in the PDR consisted mainly of 70 
cm long thermocoax cables on the sample holder. Also, 
micro- wave filters from Mini-circuits (BLP 1.9) were used 
at room-temperature. In the Leiden setup, besides the 1 m 
long thermocoax on the sample holder, additional powder 
filters (provided by Leiden Cryogenics) were present. Sim- 
ilar measurement results were achieved with both setups, 
although the effective temperature of the sample, deduced 
from the zero bias resistance, was somewhat lower with the 
Leiden setup: 45-60 mK as compared to 80-100 mK with 
the Nanoway setup (see e.g. Ref. [2"8]V 

The measurement set-up used in this work is shown 
in Fig. [3l Most of the complexity comes from the require- 
ments of the power gain measurement, which will be pre- 
sented in Sec. 4.2.4. The BOT base is DC current biased 
by a large resistor Rb, which is located at room tempera- 
ture. The size of the resistor was typically 1 — 10 Gfl. For 
differential measurements, an AC-signal is fed through a 
coupling capacitor Cc, in order to circumvent the large 
biasing resistor. Voltages are measured with low noise LI- 
75A voltage amplifiers and their output is fed into EG 
& G Instruments 7260 DSP lock-in amplifiers. An ad- 
ditional, surface- mount resistor, Rcc, is located on the 
sample holder, a few cm from the chip. Therefore, this re- 
sistor is initially also at base temperature. The role of the 
resistor is to act as load to the BOT, in order to allow a 
clear measurement of power gain. The resistor was also in 
used in noise measurements [2 7) to convert current noise 
to voltage noise. 



4 Experiment 

4.1 Observation of the Bloch nose 

Our experiments correspond to the voltage biased con- 
figurations of Fig. [TJ The voltage is gradually lifted to- 




Fig. 3. Measurement scheme for power gain, current gain and 
input impedance measurements. The BOT circuit is bounded 
by the red dashed box. The AC signal is capacitively coupled 
through the capacitor Cc- The signal from the pre-amplifiers 
LI-75A and DL 1211 are fed into lock-in amplifiers. 



wards the top of the energy band where it can overcome 
the first barrier and start to oscillate. The TV charac- 
teristics thus first shows an increase of the junction volt- 
age with zero current (when neglecting any quasiparticle 
leakage) and then at the threshold voltage where V = 
max {dE n (q)/dq} a current appears but the voltage over 
the junction goes down as the time averaged voltage tends 
to zero. When the system tunnels to a higher band, the 
Bloch oscillation is interrupted and a voltage drop is once 
again formed over the junction. Thus, for an increasing 
bias current, the voltage starts to rise again. 

The actual measurement can be done with respect to 
the external voltage V x (2-point measurement) or with re- 
spect to the real junction voltage V: either by a 4-point 
measurement |29l30j or one can simply subtract the known 
voltage Rc'I from the total (with the assumption that the 
resistor Rq is ohmic). According to the earlier transforma- 
tion between the serial and parallel configurations, we can 
compare the TV curves from our experiment with a series 
resistor with the theoretical treatment in Ref. [33] , which 
is for the parallel case by plotting I x = V x /Rc against V 
— the voltage over the J J. 

In the experiments on the Bloch Oscillating Transistor, 
we have also observed the Bloch nose and thus find the 
characteristic sign of Bloch oscillation. The IV curve of a 
Josephson junction in a resistive environment for sample 1 
(see parameters in Table 1) is shown in Fig. The curves 
are a result of a 2-point measurement where both the total 
current / and the transformation I x — V x /Rc are plot- 
ted against the total voltage V x and the junction voltage 
V = V x — IRc, respectively. In this way, one can truly see 
the Bloch nose in the form predicted in Refs. [5114) . The 
Bloch nose disappears when the ratio Ej /Eq goes down, 
a manifestation of the fact that Zener tunneling sets in 
already at lower currents. Theoretically, the voltage for 
the onset of back bending is given by Vb = 0.25e/C, when 
Ej I Ec <C 1 and the dissipation is due to quasiparticles in 
the J J [T3] • Without quasiparticle tunneling the blockade 
would become much larger, and for Ej <C Ec it would 
equal e/C. In our experiment, the observed VJ, is 23 flV, 
which is only little more than half of theoretical value: 
V b = 0.25e/C = 40 fiV. 



R. Lindell, L. Korhonen, A. Puska and P. Hakonen: Incoherent Cooper pair tunnelling and energy band dynamics 7 



v x OiV) 

20 40 60 80 

1 1 ■ 1 ■ 1 ■ 




vbiV) 

Fig. 4. IV curve of sample #1 from a 2-point measurement. 
The lowest curve shows the total current I vs. V x (for Ej/Ec = 
0.3) and the following curves are transformed to parallel con- 
figuration I x = V x /Rc vs. V = V X — IRc with Ej/E c = 0.3, 
0.24 and 0.03 from top to bottom. 

The first report on the Coulomb blockade of Cooper 
pairs was made by Haviland et al. in Ref. [3] . Bloch oscil- 
lations and the role of Zener tunneling was investigated by 
Kuzmin et al. |31|32j in experiments with a single Joseph- 
son junction in an environment of a chromium resistor. 
The characteristic Bloch nose feature, or the regions of 
negative differential resistance, was first measured for 1- 
and 2-dimensional SQUID arrays [33 34). In later exper- 
iments using SQUID arrays as a tunable, but very non- 
linear, high impedance environment, the Bloch nose could 
for the first time be clearly observed in a single Joseph- 
son junction [29 30 . Although, there the authors plot / 
vs. V and do not make the transformation I x = V x /Rc as 
done here. It should be noted, that our observation of the 
Bloch nose involves a real, linear resistive environment as 
the electromagnetic environment. 



i-i 20 




V c (mV) 



Fig. 5. IV curve for sample 1 at Ib — 40 pA (black curve) 
and 60 pA (red curve) and the DC-current gain (blue curve) 
calculated by direct subtraction of the two IV curves. The extra 
Rcc ~ 100 kfi was located at collector. T = 34 mK. 




V r (mV) 



Fig. 6. Simulated IV curve with the same device parameters 
as in the experiment presented in Figure [5] Base current 100 
pA (black curve) and 110 pA (red curve) and the DC-current 
gain (blue curve) calculated by direct subtraction of the two 
IV curves. 



4.2 Device characteristics 

4.2.1 DC Current-voltage characteristics 

We will now describe the results of experiments of sam- 
ple #1 ( parameters are given in Table [TJ whereas for the 
definition of symbols, see Fig. [3J The maximum gain was 
achieved for the maximum Ej/Ec ratio (0.3 for this sam- 
ple). Hence, if not otherwise mentioned, all the following 
presentations of the results are given for this value. Ac- 
cording to the theories on BOT operation outlined earlier, 
the gain should increase with Ej/Ec- However, the theo- 
ries assume that Ej <C Eq so that it suffices to consider 
only the dynamics between the first two bands. In the 
BOT measurements described in Ref. [8] the maximum 
gain was achieved for Ej/Ec = 3.4, which also means 
that more than the two first bands may contribute to the 
dynamics. 



In Fig. [5] a typical /U-curve of the BOT is shown for 
two values of the base current and also the corresponding 
current gain is included in the graph. The BOT has cur- 
rent amplification for both positive and negative values of 
Vc- However, as can be seen from the figure, the ampli- 
fication is largest for negative Vc- This we consider the 
" normal operating mode" of the BOT and the gain region 
for positive Vc we call the " inverted mode" . Simulated IV 
curves and current gain are shown in Fig. [6] for the same 
device parameters as in the experiment. A slightly larger 
base current than in the experiment, 100 pA as compared 
to 40 pA, was used to find a matching current gain at 
the normal operating regime. The two sets of IVs' look 
quite similar. However, the simulated one shows sharper 
features due to both a lower effective temperature and a 
higher base current, which leads to hysteretic behavior for 
the inverted operating region. Raising the temperature or 
lowering the base current in the simulation leads to a de- 
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creasing current gain. A prominent discrepancy between 
the experimental and simulated IV's is the location of 
the Bloch "back bending" and the maximum gain. In the 
experiment, this region is about 100 /iV further than in 
the simulation. This is an indication of the fact that the 
dynamics may actually involve higher bands. More sim- 
ulations for different BOT parameters can be found in 
Ref. [II]. 

4.2.2 Current gain 

Next, we consider the Ie—Ib relation for a fixed bias volt- 
age Vc and as a function of the ratio Ej/Ec- In Figs. [7] 
and [5] the Ie — Ib curves are shown for the normal and in- 
verted operating regions. In these measurements, the extra 
collector resistance Rcc was absent, hence we have larger 
operating currents. In the simulated Ie — Ib curves for 
the normal region, the gain was only about half of that 
which was observed experimentally, and the gain region is 
also much wider. However, by increasing the ratio Ej/Ec 
in the simulation to values slightly above those obtained 
from the experiment (the simulated Ej/Ec = 0.4 is in- 
cluded in Fig. [7]) , the simulated gains become closer to the 
experimental. 

The simulation for the inverted mode showed diverg- 
ing gains due to hysteretic behavior that leads to very 
low dynamic region (see also Fig. [101). This approach to 
hysteretic and divergent state could also be observed ex- 
perimentally, as noticed in Fig. [8l 




Fig. 7. Ie-Ib curves for the normal region for different Ej/Ec 
ratios. The experimental curves (solid lines) go through the 
ratios 0.03, 0.07, 0.12, 0.17, 0.21, 0.25, 0.28, 0.29 and 0.30 from 
top to bottom. The simulated curves (dotted lines) go through 
the ratios 0.25, 0.28, 0.30 and 0.40 from top to bottom. 



We measured the differential gain as a function of Vc 
with lock-in amplifiers. The results for the normal oper- 
ating region are presented in Fig. [5] for the case when the 
extra collector resistance Rcc was present. The gains are 
then a factor of 3 larger than in the above measurement 
without Rcc- This is also in close agreement to simulated 
gains, which will be discussed further in section 4.2.3. 



< 

5 




Fig. 8. Ie-Ib curves for the region of inverted operation for 
different Ej/Ec ratios. The experimental curves (solid lines) 
go through the ratios 0.03, 0.07, 0.12, 0.17, 0.21, 0.25, 0.28, 0.29 
and 0.30 from bottom to up. The simulated curves (dotted 
lines) go through the ratios 0.17, 0.21, 0.25, 0.29, 0.30 from 
bottom to top. 
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Fig. 9. The emitter current gain at the normal region for 
Ej/Ec = 0.3 and the base currents 0-100 pA, with 10 pA 
increments from right to left. The measurements were done 
using the lock-in technique of Fig. [3] Here, the 100 kS7 resistor 
at collector was present, thus leading to smaller base currents 
for optimal operation. 



The dynamic region (DR) in terms of the base cur- 
rent Ib can be inferred from the gain curves in Fig. [9] by 
using the relation: AIb = AVc/(Z out (3c)- The resulting 
DR versus flc plot is shown in Fig. QTJ The large spread in 
the datapoints are mostly due to the errors in determin- 
ing Z out by numerical differentiation of IV curves. From 
the plot we see that the DR falls as cx (3q (in fact, the 
exponent of the regression is —0.89 ± 0.12). As already 
mentioned, the simulated gains for the normal region had 
a much wider DR, though with less gain than in the ex- 
periment. 
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Fig. 10. Dynamic region (width at half maximum, in Ib- 
space) as function of j3c for the normal region for the base 
current I B = 10, 20, 50, 70 and 100 pA (from left to right). 
The solid line is a linear regression with slope —0.89 ± 0.12. 



The analytic theory in Ref. pTj predicts that the gain 
should rise exponentially with the square of the ratio Ej / Eq 



(3 C = 2 exp 



8fi \E C 



(29) 



According to our simple model presented in Sec. 2.4 
the gain dependence is non-monotonic. This could, how- 
ever, be a failure of the model as we are in a region were 
its validity is not guaranteed by the perturbation theory. 
For Ej/Ec = 0.1 the gain starts to rise again. From the 
full time-dependent simulations, the dependence is also 
exponential but less than the normal mode gains from 
experiment (Fig. ITT)) . When including heating effects in 
both the simulation and analytical model, decreased the 
gain and made the fit worse. 



4.2.3 Input impedance 

The applications of the BOT are largely determined by its 
input and output impedances, power gain and the noise 
temperature (for measurements on the noise temperature 
see [27]). We have studied the input impedance 



dVB 



(30) 



at different levels of the current gain. In Fig. [T2] the de- 
pendence is shown for the case where Ib = 60 pA. The 
input impedance and (3c were measured with lock-in am- 
plifiers with excitation frequency of 17.5 Hz. The measure- 
ment results had to be corrected for the fact that part of 
the AC-current leaks through the stray capacitance, which 
mostly originates from the thermocoax. The corrected in- 
put impedance is then 



zZ c 




0.15 0.20 0.25 

E /E 

j c 



0.30 



Fig. 11. Current gain /3c as a function of Ej/Ec- The base 
current Is was 100 pA both in the experiment and simulation. 
Experiment normal operating mode (□), simulation (dashed 
line), the analytical model presented in this paper (•) and the 
analytical approximation from Ref. 11 (thin line). The tem- 
perature was assumed to be the same as base temperature in 
.the experiment, T = 34 mK. 

where z = dv B /di B is the measured AC-differential impedance 
at the base and Zc = 1/(wCb). With the assumption that 
Cb — 600 pF, the AC and DC gains became equivalent. 

The theoretical input impedance is easily calculated 
for the "black box model" of Sec. 2.4, assuming lumped 
circuit elements and a known, constant current gain fic- 
Using the notation in Fig. [2] we can write for the island 
voltage Vi 

Vj=V c - R'lc, (32) 

where R' = Rc + Rcc, which includes the extra collector 
resistor. On the other hand we also have 

Vi = V b -IbRn. (33) 

As we know that Iq = —PcIb we get from these two 
equations 



AV B = V c + R'PcIb+IbR 
and the input impedance is then 
dV B 



N: 



dl f 



= R N + R'(3 C - 



(34) 



(35) 



(31) 



As Pc is in our definition > 0, the input impedance be- 
comes positive. The input impedance observed in the ex- 
periment (Fig. I12p shows some deviation from the simple 
model. At the point of maximum gain the impedance is 
1.4 the value of the model, but as one moves to either side 
of the maximum, the behavior is more complicated and 
also non-symmetric. The simulated curve reveals that the 
behavior could be anticipated from the dynamical model, 
which gives approximately the same 1.4 deviation from the 
simple model as observed in the experiment. The asym- 
metric feature shows that the simple model, that Zi„ is 
proportional to /3c, is followed quite closely for the part of 
the current gain where voltage is higher than the optimum 
operating point. 
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Pc 

Fig. 12. Z in /(f3 c R') as function of (3 C for I B = 60 pA for 
experiment (blue curve) and simulated case (red curve). The 
lower part of the curves is the behavior left of the maximum 
gain and the upper part is for the right side. Inset: Zi„ and f3c 
as function of Vc for experiment (thick lines) and simulation 
(thin lines). 

4.2.4 Power gain 

An important measure of the device performance is the 
power gain 

V = Pout /Pin, (36) 

where Pj„ = i 2 B Zi n and P ou t = i c Z out , and the currents 
are rms values. The simple black-box model for the input 
impedance, where Zi n ~ (3c Re, gives a theoretical power 
gain of 

V=^Pc (37) 

The setup for the power gain measurement is shown 
in Fig. [31 The AC-signal was capacitively coupled to by- 
pass the large bias resistor Rb- We measured the power 
delivered by the BOT to a 100 kSl resistor (Rcc), which 
acted as a load at the collector. The measured power at 
the output is given by P ou t = v cc/-^cc, where vec is 
the measured AC- voltage over the load. The measured 
power gain in this case becomes r\ — {v cc / 'i 2 B ) / {Rcc Zi n ) ■ 
The simple black box model gives a power gain of 77 = 
{Rcc I R')l3c- Taking into account the measured input 
impedance in Fig. [12] we find that, for Ib = 60 pA, 77 ~ 
{Rcc /l-4i£')/3c — 6, which agrees with the independent 
measurement in Fig. 1131 

The power gains for Ib = — 100 pA at the normal 
operating point are shown in Fig. [13] The largest measured 
power gain was around 35 for Ib = 100 pA. The output 
impedance of the device itself at the operating point was in 
the range —50 kSl to —100 kO. The gain might be expected 
to grow rapidly when Rcc approaches \Z out \. However, a 
positive load impedance with similar absolute magnitude 
as the negative output impedance may lead to oscillatory 
behavior. 

A simulated power gain for Ib — 60 pA is shown in the 
inset. The agreement looks fairly good but when the base 



current is increased in the simulation the effect on the gain 
is quite small compared to the large increase observed in 
the experiment. The simulation was also quite sensitive 
to the capacitance Ccc (see Fig. [3j , here we settled for a 
ratio Cb I Ccc = 50, which gave a good balance between 
achieving a more realistic model and simulation stability. 
The value of Cb used in the simulation was 1 pF, which 
means that the base voltage drops less than 0.1 % when 
a quasiparticle tunnels to the island. 




-400 -380 -360 

Fig. 13. Power gain with 100 kQ. load for I B =0, 20, 60, 80 
and 100 pA from right to left and Inset: Power gain for Ib = 60 
pA (•) and simulated curve (red line) 



5 Conclusions 

We have here briefly reviewed the physical principles and 
computational methods that have been applied to ana- 
lyze the Bloch Oscillating Transistor. The circuit has been 
shown to produce a variety of interesting physical phenom- 
ena that also have a potential for practical applications. 
The measurements on the BOT have shown that the de- 
vice works according to the physical principles discussed. 
The Bloch nose as a manifestation of the competition be- 
tween coherent current via Bloch oscillations and Coulomb 
blockade of Cooper pairs was observed in a controlled re- 
sistive environment. The smallness of the observed block- 
ade would indicate that the environmental dissipation is 
dominated by quasiparticles. 

The optimal BOT parameters are quite difficult to 
specify, mainly, because the used models are approxima- 
tions that hold in a certain regime for the parameters. We 
seek to optimize Rn, Rjj, R, Ej/Ec and the capacitance 
C b , and at the same time keeping in mind what the opti- 
mum input and output impedances should be for a specific 
application. The analytic model of Sec 2.4 gives a closed 
formula for the current gain, and can thus used for opti- 
mization. In all the models, the gain grows with Ej/Ec, 
but simultaneously, we have the requirement Ej/Ec <C 1 
for the theories to be valid. However, nothing stops us to 
go to somewhat larger Ej/Ec ratios in the simulation, 
although, the bands then increasingly deviate from the 
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parabolic approximation. The simulation shows that the 
gains also increase with Ej/Eq, but naturally, with the 
expense of the dynamic region. The simple analytic model 
does not depend on Rn at all. However, in the analyt- 
ical model of Ref. [TT] the device performance improves 
as Rn is reduced. A too transparent base junction would, 
however, be undesirable for many reasons: the increased 
transparency could increase charge fluctuations and de- 
stroy the Bloch oscillations in the Joscphson junction. In 
fact, in our samples the NIS junction was quite resistive. 
Also, small Rn would, in practice, mean larger capaci- 
tance and smaller blockade. The model indicates that a 
larger collector resistance Rc is always desirable. This is 
probably a general conclusion, bearing in mind that the 
input impedance increases linearly with Rc- In practice, 
fabricating chromium on-chip resistors larger than 0.5 Mf2 
has proven difficult due to nonlinearities. 

The observed device properties of the BOT can be 
qualitatively reproduced by simulations and the simple 
theoretical considerations presented here. There are, how- 
ever, still quite many discrepancies between simulations 
and experiment. Often the values of the simulation for 
gains and impedances were a factor 2 from the experimen- 
tal case. The simulations do show that the main operating 
principle of the BOT can be understood by the processes 
outlined. There are still some places for improvement, 
e.g., by taking into account the true phase fluctuations 
caused by tunneling of quasiparticles to the island. But, 
this would mean a great increase in the complexity of the 
simulation, introducing new numerical challenges. Also, 
taking into account the true band structure of the device 
might improve the agreement, especially for Ej > Eq- 
While awaiting a more complete theoretical treatment, 
the simulations can be used for qualitative modelling of 
the BOT. The device could still be used as an on-chip 
amplifier and detector for mesoscopic experiments, where 
a good current gain is needed at medium level impedances, 
but, the exact device properties can be inferred only from 
experiment. 
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Appendix 

Transconductance 

The transconductance can be calculated in the simple 
black box model by g = dlc/dVB = —f3c(dlB I 'dVs) = 
-Pc/(Rcf3c) = -l/Rc- For sample #1, 1/R C = 5.3 
and the maximum observed transconductance was 5.0±0.5 

us. 



For sample #2, which had a larger Rc = 368 kfl, the 
transconductance was 9.1 fiS, which is about 3.5 larger 
than expected from the simple model above (see Fig. [T4|) . 
The IVs cross at the dips in the TV-curves, thus giving 
rise to a sign change in transconductance, and also for 
the current gain. Thus, the input impedance Zi n — (3c I g 
stays positive in all measurements. 




F c (mV) 

Fig. 14. IV and transconductance curve for BOT #2. Vb = 
(black curve) and Vb = 17.9 /iV (red curve). The dotted lines 
show the points where the IVs cross and the transconductance 
changes sign. 

The IV curve for sample #2 was intrinsically hys- 
teretic in behavior when applying a current bias to the 
base via an additional resistor at room-temperature. The 
hysteresis, however, disappeared when the biasing resistor 
was between 100 kf2 and 1 MO, which could indicate that 
the device has a region of negative input impedance in this 
range. Simulations have shown that the input impedance 
is always positive, but that the base Ib — Vb curve can be 
hysteretic and, therefore, the observed hysteresis is a con- 
sequence of the device switching between two operating 
regimes, which have the same base current but different 
impedance levels. 
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